Spatial scales of COVID-19 transmission in Mexico (preprint)

During outbreaks of emerging infectious diseases, internationally connected cities often experience large and early outbreaks, while rural regions follow after some delay. This hierarchical structure of disease spread is influenced primarily by the multiscale structure of human mobility. However, during the COVID-19 epidemic, public health responses typically did not take into consideration the explicit spatial structure of human mobility when designing non-pharmaceutical interventions (NPIs). NPIs were applied primarily at national or regional scales. Here we use weekly anonymized and aggregated human mobility data and spatially highly resolved data on COVID-19 cases, deaths and hospitalizations at the municipality level in Mexico to investigate how behavioural changes in response to the pandemic have altered the spatial scales of transmission and interventions during its first wave (March - June 2020). We find that the epidemic dynamics in Mexico were initially driven by SARS-CoV-2 exports from Mexico State and Mexico City, where early outbreaks occurred. The mobility network shifted after the implementation of interventions in late March 2020, and the mobility network communities became more disjointed while epidemics in these communities became increasingly synchronised. Our results provide actionable and dynamic insights into how to use network science and epidemiological modelling to inform the spatial scale at which interventions are most impactful in mitigating the spread of COVID-19 and infectious diseases in general.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico (preprint)

Up to November 2021, over 200 different SARS-CoV-2 lineages circulated in Mexico. To investigate lineage replacement dynamics, we applied a phylodynamic approach to explore the evolutionary trajectories of five dominant lineages that circulated during the first year of the local epidemic. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in the country. Lineages B.1.1.222 and B.1.1.519 showed comparable dynamics, represented by clades likely originating in Mexico and persisting for over a year. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. We further explored viral movements across the country, applied within the largest clades identified (belonging to lineage B.1.617.2). Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico (preprint)

Up to November 2021, over 200 different SARS-CoV-2 lineages circulated in Mexico. To investigate lineage replacement dynamics, we applied a phylodynamic approach to explore the evolutionary trajectories of five dominant lineages that circulated during the first year of the local epidemic. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in the country. Lineages B.1.1.222 and B.1.1.519 showed comparable dynamics, represented by clades likely originating in Mexico and persisting for over a year. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. We further explored viral movements across the country, applied within the largest clades identified (belonging to lineage B.1.617.2). Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Emergence and widespread circulation of a recombinant SARS-CoV-2 lineage in North America (preprint)

Genetic recombination is an important driving force of coronavirus evolution. While some degree of virus recombination has been reported during the COVID-19 pandemic, previously detected recombinant lineages of SARS-CoV-2 have shown limited circulation and been observed only in restricted areas. Prompted by reports of unusual genetic similarities among several Pango lineages detected mainly in North and Central America, we present a detailed phylogenetic analysis of four SARS-CoV-2 lineages (B.1.627, B.1.628, B.1.631 and B.1.634) in order to investigate the possibility of virus recombination among them. Two of these lineages, B.1.628 and B.1.631, are split into two distinct clusters (here named major and minor). Our phylogenetic and recombination analyses of these lineages find well-supported phylogenetic differences between the Orf1ab region and the rest of the genome (S protein and remaining reading frames). The lineages also contain several deletions in the NSP6, Orf3a and S proteins that can augment reconstruction of reliable evolutionary histories. By reconciling the deletions and phylogenetic data, we conclude that the B.1.628 major cluster originated from a recombination event between a B.1.631 major virus and a lineage B.1.634 virus. This scenario inferred from genetic data is supported by the spatial and temporal distribution of the three lineages, which all co-circulated in the USA and Mexico during 2021, suggesting this region is where the recombination event took place. We therefore support the designation of the B.1.628 major cluster as recombinant lineage XB in the Pango nomenclature. The widespread circulation of lineage XB across multiple countries over a longer timespan than the previously designated recombinant XA lineage raises important questions regarding the role and potential effects of recombination on the evolution of SARS-CoV-2 during the ongoing COVID-19 pandemic.

Identification of site-specific evolutionary trajectories shared across human betacoronaviruses (preprint)

Comparison of evolution among related viruses can provide insights into shared adaptive processes, for example following host switching to a mutual host species. Whilst phylogenetic methods can help identify mutations that may be important for evolutionary processes such as adaptation to a new host, these can be enhanced by positioning candidate mutations to known functional sites on protein structures. Over the past two decades, three zoonotic betacoronaviruses have significantly impacted human public health: SARS-CoV-1, MERS-CoV and SARS-CoV-2, whilst two other betacoronaviruses, HKU1 and OC43, have circulated endemically in the human population for over 100 years. In this study, we use a comparative approach to prospectively search for potentially evolutionarily-relevant mutations within the Orf1ab and S genes across betacoronavirus species that have demonstrated sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1 and SARS-CoV-2). We used a combination of molecular evolution methods to identify 30 sites that display evidence of homoplasy and/or stepwise evolution, that may be suggestive of adaptation across emerging and endemic betacoronaviruses. Of these, seven sites also display evidence of being selectively relevant. Drawing upon known protein structure data, we find that four of the identified mutations [18121 (exonuclease/27), 21623 (spike/21), 21635 (spike/25) and 23948 (spike/796), in SARS-CoV-2 genome coordinates] are proximal to regions of known functionality. Our results provide a molecular-level context for common evolutionary pathways that betacoronaviruses may undergo during adaptation to the human host.

Genomic epidemiology of SARS-CoV-2 transmission lineages in Ecuador (preprint)

Characterisation of SARS-CoV-2 genetic diversity through space and time can reveal trends in virus importation and domestic circulation, and permit the exploration of questions regarding the early transmission dynamics. Here we present a detailed description of SARS-CoV-2 genomic epidemiology in Ecuador, one of the hardest hit countries during the early stages of the COVID-19 pandemic. We generate and analyse 160 whole genome sequences sampled from all provinces of Ecuador in 2020. Molecular clock and phylgeographic analysis of these sequences in the context of global SARS-CoV-2 diversity enable us to identify and characterise individual transmission lineages within Ecuador, explore their spatiotemporal distributions, and consider their introduction and domestic circulation. Our results reveal a pattern of multiple international importations across the country, with apparent differences between key provinces. Transmission lineages were mostly introduced before the implementation of non-pharmaceutical interventions (NPIs), with differential degrees of persistence and national dissemination.

Genomic analysis of early SARS-CoV-2 strains introduced in Mexico (preprint)

The COVID-19 pandemic has affected most countries in the world. Studying the evolution and transmission patterns in different countries is crucial to implement effective strategies for disease control and prevention. In this work, we present the full genome sequence for 17 SARS-CoV-2 isolates corresponding to the earliest sampled cases in Mexico. Global and local phylogenomics, coupled with mutational analysis, consistently revealed that these viral sequences are distributed within 2 known lineages, the SARS-CoV-2 lineage A/G, containing mostly sequences from North America, and the lineage B/S containing mainly sequences from Europe. Based on the exposure history of the cases and on the phylogenomic analysis, we characterized fourteen independent introduction events. Additionally, three cases with no travel history were identified. We found evidence that two of these cases represent local transmission cases occurring in Mexico during mid-March 2020, denoting the earliest events described in the country. Within this Mexican cluster, we also identified an H49Y amino acid change in the spike protein. This mutation is a homoplasy occurring independently through time and space, and may function as a molecular marker to follow on any further spread of these viral variants throughout the country. Our results depict the general picture of the SARS-CoV-2 variants introduced at the beginning of the outbreak in Mexico, setting the foundation for future surveillance efforts. This work is the result of the collaboration of five institutions into one research consortium: three public health institutes and two universities. From the beginning of this work, it was agreed that the experimental leader of each institution would share the first authorship. Those were the criteria followed to assign first co-first authorship in this manuscript. The order of the other authors was randomly assigned. IMPORTANCEUnderstanding the introduction, spread and establishment of SARS-CoV-2 within distinct human populations is crucial to implement effective control strategies as well as the evolution of the pandemics. In this work, we describe that the initial virus strains introduced in Mexico came from Europe and the United States and the virus was circulating locally in the country as early as mid-March. We also found evidence for early local transmission of strains having the mutation H49Y in the Spike protein, that could be further used as a molecular marker to follow viral spread within the country and the region.